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Abstract 

The diffractive photoproduction of J/ifi mesons is measured with the HI detector at the ep 
collider HERA using an integrated luminosity of 78 pb _1 . The differential cross section 
d<r(7p — > J/ipY)/dt is studied in the range 2 < \t\ < 30 GeV 2 , where t is the square of 
the four-momentum transferred at the proton vertex. The cross section is also presented as 
a function of the photon-proton centre-of-mass energy W 1V in three t intervals, spanning 
the range 50 < W 1P < 200 GeV. A fast rise of the cross section with W JP is observed 
for each t range and the slope for the effective linear Pomeron trajectory is measured to be 
a' = -0.0135 ±0.0074 (stat.) ±0.0051 (syst.) GeV" 2 . The measurements are compared 
with perturbative QCD models based on BFKL and DGLAP evolution. The data are found 
to be compatible with s-channel helicity conservation. 
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1 Introduction 



The diffractive photoproduction of J / ip mesons with large negative momentum transfer squared 
t at the proton vertex is a powerful means to probe the parton dynamics of the diffractive ex- 
change. The variable t provides a relevant scale to investigate the application of perturbative 
QCD (pQCD). The diffractive photoproduction of vector mesons can be modelled in the proton 
rest frame, where the photon fluctuates into a qq pair at a long distance from the proton target. 
The colour singlet exchange between the qq fluctuation and the proton is realised in lowest or- 
der QCD by the exchange of a pair of gluons with opposite colour. In the leading logarithmic 
(LL) approximation, this process is described by the effective exchange of a gluonic ladder. At 
sufficiently low values of Bjorken x (i.e. large values of the centre-of-mass energy W 1P ), the 
gluon ladder is expected to include contributions from BFKL evolution [1|, as well as from 
standard DGLAP evolution Q. Compared with other channels which have been used to search 
for BFKL evolution 001 El El Eld, the measurement of diffractive J/ip production at large 
\t\ provides an experimentally clean signature in which the accurate measurement of the J/il> 
four-momentum allows the kinematic dependences of the process to be determined precisely. 

In this paper, an analysis of the diffractive photoproduction process 'jp — > J / ipY is presented, 
extending into the hitherto unexplored region of large \t\ (2 < \t\ < 30 GeV 2 ). Here, the system 
Y represents either an elastically scattered proton or a dissociated proton system. For the range 
of \t\ studied in this analysis, the contribution from elastic J/ip production may be neglected 
due to its steep \t\ dependence 0. The cross section is measured differentially as a function of 
\t\ and as a function of the photon-proton centre-of-mass energy W 1P in different regions of \t\, 
using the J/ip decay into two muons. To obtain information about the helicity structure of the 
interaction, the spin density matrix elements are extracted. 



2 Perturbative QCD Models 

Perturbative QCD models for the photoproduction of J/ip mesons have been developed in the 
leading logarithmic approximation using either BFKL ifTDlfTTllT^l or DGLAP lfT51 evolution. 
In the BFKL LL model the cross section depends linearly on the parton distribution of the 
proton and the gluon ladder couples to a single parton (dominantly a gluon) within the proton. 
The BFKL amplitude is expanded in terms of log(x/j, W^ p /Wq), where Xh is the fraction of the 
proton momentum carried by the parton struck by the diffractive exchange. The scale parameter 
Wq is chosen to be half the vector meson mass My. The value of a s is fixed in the model 
to a value consistent with that extracted from a fit lTl"2l to proton dissociative p, <p and J/ip 
photoproduction data at HERA lfT4l . The BFKL LL model predicts an approximate power-law 
behaviour for the t dependence of the form der/dt oc \t\~ n , where n is a function of \t\. For 
the kinematic range studied here, n increases from around 3 to 4 with increasing |t| and the 
approximation to a power-law improves as \t\ increases. The calculation predicts a fast rise of 
the cross section a ~ with 5 ~ 1.4, which has little or no dependence on the value of t. In 
a recent paper lfi~5l . the LL calculations have been extended to incorporate the effects of higher 
conformal spin lfT*6l . Although the full next-to-leading order terms of the BFKL amplitude have 
yet to be calculated for non-zero t, an estimate of the non-leading (NL) corrections was obtained 
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using kinematic constraints. In the DGLAP LL model, the cross section depends on the squared 
gluon distribution of the proton. The model predicts a non-exponential t dependence and a steep 
energy dependence which flattens as \t\ approaches M v due to the limited phase space available 
for evolution. 

In the pQCD models EH QUEH [HUH), a non-relativistic approximation [17] for the J/i/j 
wavefunction is used in which the longitudinal momentum of the vector meson is shared equally 
between the quark and the anti-quark. In this approximation, the vector meson retains the 
helicity of the photon such that s-channel helicity conservation (SCHC) is satisfied lfl~8l . 



3 Data Analysis 

The data presented here were recorded in the years 1996 to 2000 and correspond to an integrated 
luminosity of 78 pb -1 . The majority of the data were collected when HERA was operated with 
positrons of energy 27.5 GeV and protons of 920 GeV. These data are combined with smaller 
data samples in which either the proton energy was 820 GeV or the positrons were replaced by 
electrons. 



3.1 The HI Detector 

A detailed description of the HI detector can be found in lfi"9l and only a short overview of 
the detector components most relevant to the present analysis is given here. The z-axis of the 
HI detector is defined along the beam direction such that positive z values correspond to the 
direction of the outgoing proton beam. 

Charged particles emerging from the interaction region are measured by the central tracking 
detector (CTD) in the pseudorapidity range —1.74 < r] < 1.74 1 . The CTD comprises two 
large cylindrical central jet drift chambers (CJC) and two z-chambers arranged concentrically 
around the beam-line within a solenoidal magnetic field of 1.15 T. The CTD also provides 
triggering information based on track segments in the r — (j> plane from the CJC and the z- 
position of the vertex from a double layer of multi-wire proportional chambers. The energies 
of final state particles are measured in the liquid argon (LAr) calorimeter, which surrounds the 
tracking chambers and covers the range —1.5 < r] < 3.4. The backward region (—4.0 < r] < 
—1.4) is covered by a lead-scintillating fibre calorimeter (SPACAL [ 20 1) with electromagnetic 
and hadronic sections. The calorimeters are surrounded by the iron return yoke of the solenoidal 
magnet. The tracks of muons which penetrate the main detector are reconstructed from streamer 
tubes placed within the iron in the range —2.5 < r] < 3.4. The luminosity is measured using 
the small angle Bremsstrahlung process (ep — > ep'j) in which the final state photon is detected 
in a calorimeter, close to the beam-pipe, at 103 m from the nominal interaction point. 

'The pseudorapidity ?/ of an object detected with polar angle 9 is defined as i] = — In tan(#/2). 
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3.2 Kinematics 



The kinematics for diffractive charmonium production ep — > eJ/ipY are described in terms of 
the ep centre-of-mass-energy squared s = (k + p) 2 , the virtuality of the photon 
Q 2 = —q 2 = —(k — k') 2 , the square of the centre-of-mass energy of the initial photon-proton 
system W 2 p = (q + p) 2 and the four-momentum transfer squared t — (j> — py) 2 . Here k (k') 
is the four-momentum of the incident (scattered) lepton and q is the four-momentum of the 
virtual photon. The four-momentum of the incident proton is denoted by p and py is the four- 
momentum of the system Y. The event elasticity is defined as z = (p ■ p^) /(p ■ q) where p^ is 
the four-momentum of the J/if; . In the proton rest frame z is equal to the fractional energy of 
the photon transferred to the vector meson. 



3.3 Event Selection 

In this analysis, the J/ip mesons are detected via their decay into two oppositely charged muons 
(branching fraction 5.88 ± 0.10% lETl ). The data were selected by a combination of triggers 
based on muon and track signatures. The selected events are required to have a vertex located 
in z within 40 cm of the nominal interaction point. Events with two tracks of opposite charge 
in the CJC, each associated with the event vertex and each with pseudo-rapidity \rj\ < 1.7 '4 and 
transverse momentum p T > 0.8 GeV are used to form J/xj) candidates. Both decay muons are 
identified in the instrumented iron or as minimum ionising particles in the LAr calorimeter. 

Photoproduction events are selected by the absence of a scattered beam lepton candidate in the 
LAr or SPACAL calorimeters. The accepted photoproduction event sample covers the range 
Q 2 < 1 GeV 2 with an average (Q 2 ) ~ 0.06 GeV 2 , as determined from Monte Carlo simula- 
tions. 

In order to select diffractive events, the analysis is restricted to the region of elasticity z > 0.95. 
For the range of t and W lp studied in this paper, the cut z > 0.95 restricts the invariant mass of 
the system Y to be in the range M Y < 30 GeV, through the relation z ~ 1 — (M Y — t) /W 2 p . 
The measurement of z is obtained from (E —p z )j/^,/ J2(E~Pz) wnere XX Pz) is calculated 
from all detected particles in the calorimeters and the CJC including the decay products of the 
J/ip. The variable W yp is reconstructed using W 2 p = J2(E ~ Vz) • %E P where E p is the energy 
of the incident proton beam. In the kinematic region studied, the variable t is well approximated 
by the negative transverse momentum squared of the vector meson, i.e. t ~ —p 2 tJ i^- 



3.4 Monte Carlo Simulation 

Monte Carlo simulations are used to correct the data for the effects of resolution, acceptance and 
efficiency losses. Samples of events from signal and background processes are passed through 
a detailed simulation of the detector response, based on the GEANT program lf22l . and through 
the same reconstruction software as was used for the data. 

The Monte Carlo generator used for the simulation of proton dissociative diffractive J/tp pro- 
duction is HITVM ll23l . which generates events according to the BFKL model described in 
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HlOHl II . The events are generated using the GRV94-HO partem density functions ll2"4l and the 
partonic system is fragmented according to the Lund string model implemented within the JET- 
SET program E51 . The generated M Y distribution in HITVM has an approximate exponential 
dependence du/dMy ~ e -o.iA/ y SCHC is assumed for the photon to vector meson transition. 

The final sample of events contains background from resonant and non-resonant sources. The 
resonant background is produced indirectly through the decay of ip(2S) mesons. This contribu- 
tion is simulated using a Monte Carlo sample of ip(2S) mesons generated using the DIFFVM 
Monte Carlo generator Eol according to the 4>(2S) t distribution and cross section ratio to J/ip 
production measured at lower values of \t\ 0. A contribution of 4% is observed with no sig- 
nificant t dependence. The main contribution to the non-resonant background is from the QED 
77 — > fifj, process, which is simulated using the LPAIR ETll Monte Carlo generator. 

3.5 Signal Extraction 

The invariant mass spectrum for all events in the range \t\ > 2 GeV 2 , 50 < W 1P < 150 GeV and 
z > 0.95 is shown in figured] The LPAIR non-resonant background is normalised to the data 
in the side-bands outside the mass regions of the J/ip and ip(2S) resonances. The number of 
signal events is determined from the number of events in the mass window of 2.9 < Af M + M - < 
3.3 GeV, after subtracting the contributions of the resonant and non-resonant backgrounds. 
The resulting number of J ftp candidate events for the total sample shown in figure [T] is 846 ± 
30 (stat.). 

3.6 Comparison of Data and Simulation 

The HITVM model gives a reasonable description of the data which is further improved through 
small adjustments to the W yp and t distributions. After these adjustments a comparison between 
the simulation and the data, before background subtraction, is given in figure El for the region 
\t\ > 2 GeV 2 , 50 < W 1P < 150 GeV, z > 0.95 and 2.9 < M M+M - < 3.3 GeV. Distributions 
are shown for the polar angle and transverse momentum of the decay muon tracks, for the 
reconstructed value of the elasticity z (where the cut on z is not applied), for W 1P , for the decay 
angular distributions cos#* and 0* (see section l4~2l and for the squared transverse momentum 
of the dimuon system p 2 t ^ +fjL _. The structure in the 0* distribution (figure Ef) is due to the 
low acceptance for one of the muons, which has a low transverse momentum in the laboratory 
frame, when the J/ip meson production and decay planes coincide (<p* ~ 0° or (ff ~ ±180°). 

3.7 Systematic Uncertainties 

The uncertainties in detector effects and in the modelling of the underlying physics processes 
contribute to the systematic uncertainties in the cross section measurements. The following 
sources of systematic error are taken into account. 



7 



• The uncertainty in the acceptance corrections is estimated by reweighting the W 1P distri- 
bution by W^p- 35 and the t distribution by t ±0 - 85 . The resulting systematic uncertainties 
on the cross section measurements range from 1% to 5%. 

• The uncertainty in the mass distribution of the proton dissociative system Y is estimated 
by reweighting the My dependence in HITVM by e ±omMY . This results in a variation of 
the cross section of about 4%, increasing up to 19% at the largest W 1P and \t\. 

• The effect of possible deviations from SCHC is estimated by modifying the simulated 
cos 9* distribution. The cross sections alter by 5% on average. 

• The uncertainty on the trigger efficiency, obtained from an independently triggered sam- 
ple of events, gives a contribution to the systematic error of 6%. 

• The uncertainty in the identification efficiency of muons is estimated by detailed compar- 
ison of the data and simulation efficiencies for an independent data sample. The resulting 
systematic uncertainty is 6%. 

• The uncertainty due to the reconstruction efficiency of the central tracker for the two 
tracks leads to an error of 4%. 

• The uncertainty in the non-resonant background subtraction is estimated by using a data 
side-band subtraction as an alternative to the Monte Carlo subtraction. A difference of 
~ 2% is found between the two methods and assigned to the systematic error. 

• The uncertainty in the subtraction of the ip(2S) background leads to an error of 2%, 
obtained by varying the normalisation and exponential t slope of the ipi^S) cross section 
in the simulation. 

• Other sources of systematic error are the uncertainty in the hadronic energy scale of the 
liquid argon calorimeter, the uncertainty in the luminosity measurement and the uncer- 
tainty in the branching fraction for the measured decay channel ll2"Tll . Each of them is 
responsible for an error of no more than 1.7%. 

The total systematic error for each data point has been obtained by adding all individual contri- 
butions in quadrature. It has a small dependence on t with an average value of 12% and increases 
from around 11% at low W 1P to 20% at high W JP . The part of the uncertainty which is uncorre- 
cted between different data points contributes 8.5% to the systematic error. The statistical error 
is larger than the systematic error in the region \t\ > 5.5 GeV 2 . 



4 Results 

4.1 Cross Sections 

The differential cross section da/dt for the process ep — > eJ/ipY is obtained from the number 
of data events in each measurement interval after corrections for backgrounds and detector 
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effects, divided by the integrated luminosity, the branching fraction and the width of the interval. 
The cross section for the photoproduction process 7J9 — > J/ipY is obtained by dividing the 
differential ep cross section by the effective photon flux H28I integrated over the W 1P and Q 2 
ranges of the measurement. QED radiative effects are estimated to be less than 1% and are 
neglected. The differential photoproduction cross section da/dt is shown in figure[3]and table[l] 
for the kinematic region 50 < W 1P < 150 GeV and z > 0.95. The data are plotted at the 
mean value in each t interval according to a parameterisation of the data. In the region |t| > 
3.45 GeV 2 , the data in figure are adequately described by a power-law dependence of the 
form A • \t\~ n where n = 3.00 ±0.08 (stat.) ±0.05 (syst.). When the power-law fit is repeated, 
each time increasing the starting value of \t\ in the fit, the value of n is found to increase 
systematically up to a value of n = 3.78 ± 0.17 (stat.) ± 0.06 (syst.) for \t\ > 10.4 GeV 2 . 
The data are incompatible with an exponential behaviour da/dt oc e bt which was found to 
give a reasonable description of the proton dissociative J/xj) cross section at lower values of \t\ 
< 5 GeV 2 ) 0. 

In figure the data are compared with the predictions from pQCD calculations in the BFKL 
leading logarithmic approximation [TT31 (solid curve), including non-leading corrections with 
fixed a s lfl~5l (dashed curve) and including non-leading corrections with running a s f\5\ (dot- 
ted curve). The t dependence and normalisation of the data are well described by the BFKL 
LL approximation when the parameters of the model are set to values consistent with those 
extracted from a fit lfT2l to various vector meson proton dissociation data at HERA covering a 
smaller \t\ range [14], i.e. the scale parameter is set to W = My/2 and a s is fixed at 0.18. The 
normalisation uncertainty due to the choice of Wo is large. For example, using Wq = My /4 
(Wo = My) leads to an increase (decrease) in the normalisation of the prediction by a factor 
of approximately two. The inclusion of NL corrections with a fixed strong coupling a s leads 
to only a small difference with respect to the LL prediction. However, with a running a s the t 
dependence becomes steeper and the prediction is unable to describe the data across the whole 
t range. The uncertainties in the choice of the scale parameter, proton parton density and other 
parameters used in the NL calculation have only a small effect on the shape of the predictions in 
comparison to the treatment of a s . The data are also compared with calculations in the DGLAP 
LL approximation ifTSll (dashed-dotted curve) in the region of validity for the model \t\ < Mj,^. 
The data are well described in shape and normalisation when the separation parameter to, which 
represents the value of t at which the prediction for proton dissociation matches the elastic cross 
section, is set to —0.60 GeV 2 . 

The ZEUS collaboration has recently published data on the diffractive production of J/ip 
mesons with proton dissociation in the range 1.2 < \t\ < 6.5 GeV 2 , 80 < W 1P < 120 GeV and 
x h = N/(W 2 p (l — z)) > 0.01 iTBl . When the present analysis is performed in this kinematic 
region, good agreement between the HI and ZEUS results is observed. 

In figure HJand tables El- 111 the cross section a 1P _j/^ Y is presented as a function of W lp for 
three ranges of t in the kinematic region z > 0.95. The data in each t range are consistent 
with a power-law dependence of the form a oc W^ p and the results of power-law fits for S are 
given in table |5J The contribution from correlated systematic errors is calculated by shifting 
the data points according to each source of uncertainty and repeating the fits. The values of the 
power 5 in each t range are similar to the results from the proton elastic process for J/tp mesons 
at low |£| measured over a similar range of W 1P ll29l . In a Regge pole model, the power-law 
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dependence can be expressed as da/dt = F(t)W^p ' 4 where F(t) is a function of t only. 
The value of a(t) at each t value is obtained from a = (5 + 4)/4 and is also shown in table 
Assuming a single effective Pomeron trajectory of the linear form a(t) = a(0) + a't, a fit to the 
three a values yields a slope of a' = -0.0135 ± 0.0074 (stat.) ± 0.0051 (syst.) GeV~ 2 with 
an intercept of a(0) = 1.167 ± 0.048 (stat.) ± 0.024(syst.). The value of the slope parameter 
a' is lower than that observed for the elastic photoproduction of J/i/j mesons at low \t\ O0"l . It 
is also significantly different from observations at low \t\ in hadron-hadron scattering, where a 
value of a' = 0.26 ± 0.02 GeV~ 2 El was obtained. 

In figure|Hthe data are compared with the BFKL theoretical predictions for the LL approxima- 
tion (solid curve) and the LL+NL prediction with fixed a s (dashed curve). The data are also 
compared with the DGLAP LL predictions (dashed-dotted curve). The BFKL LL contribution 
gives a reasonable description of the energy dependence, except for the lowest \t\ range where 
it is steeper than the data. The BFKL LL+NL prediction with fixed a s is similar to that of 
the BFKL LL prediction. The DGLAP LL model, which is valid in the range \t\ < Mj,^, 
describes the energy dependence in the lowest \t\ range, 2 < \t\ < 5 GeV 2 . In the region 
5 < \t\ < 10 GeV 2 , where \t\ approaches Mj,,, the description becomes worse. 



4.2 Spin Density Matrix Elements 

The polar (#*) and azimuthal (0*) decay angular distributions are measured in the rest frame 
of the J/xj) with the quantisation axis taken as the direction of the meson in the photon-proton 
centre-of-mass frame (helicity frame). The normalised two-dimensional angular distribution 
for the decay of the J/tp meson to fermions is written in terms of spin density matrix elements 
r^, r?ii and Re{r%} 021 as 

VO-^-lWr (1) 



adcosfl*d0* 4vr \2 y uu/ 2 

+v / 2Re{r^} sin 29* cos 0* + sin 2 9* cos 

The one-dimensional distributions are obtained by integrating over cos 9* or <p* and give d pse* (X 
l + r oo + (l — 3?"qq) cos 2 9* and ^ oc cos 20*. Under the assumption of s-channel helic- 

ity conservation (SCHC), the J/tp meson in photoproduction is expected to be fully transversely 
polarised and the matrix elements Tqq, r®ti and Rejr^} are zero. 

The spin density matrix elements are extracted by a two-dimensional log likelihood fit of the 
data to equation ((TJ. The normalised single differential distributions in cos 9* and <fi* are shown 
in figure |5] for three ranges of t. The dashed curve on the figure shows the expectation from 
SCHC and the solid curves show the results of the two-dimensional fit. The values of the three 
extracted matrix elements are shown in figure |6] and table |6] as a function of \t\. Measurements 
from the ZEUS collaboration of the spin density matrix elements for the photoproduction of p° 
and J/ij) mesons lfT4l are also shown in the figure. In contrast to the p° meson, the measured 
spin density matrix elements of the J/ip meson are all compatible with zero, within experi- 
mental errors, and are thus compatible with SCHC. The J/ip results are therefore consistent 
with the longitudinal momentum of the photon being shared symmetrically between the heavy 
quarks. Hence, the approximations made in the pQCD models H101llllfr21ll3lfT5l for the J/ip 
wavefunction are satisfactory for the present data. 
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5 Summary 



The differential cross section da / dt for the diffractive photoproduction of J/ tp mesons has been 
measured as a function of the momentum transfer squared t from \t\ = 2 GeV 2 up to values 
as large as \t\ = 30 GeV 2 in the kinematic region z > 0.95 and 50 < W yp < 150 GeV. The 
data are well described in this region by pQCD calculations lfl"5l using the leading logarithmic 
BFKL equation with parameters consistent with a fit to vector meson proton dissociation data 
at HERA ifm. The addition of non-leading corrections preserves the description of the data if 
the strong coupling a s is held fixed. The data in the region \t\ < Mj,^ are well described by a 
model [TT~3l based on DGLAP evolution. 

The cross section has also been measured as a function of W 1P in three t intervals. The energy 
dependence shows a similar steep rise to that observed for elastic J/ip production at low \t\ ll29l 
and the rise persists to the largest \t\ values studied. The energy dependence is reasonably 
described by the BFKL model with the chosen parameters, except for the lowest \t\ range (|t| < 
5 GeV 2 ). The DGLAP model describes the energy dependence in the range |t| < 5 GeV 2 . 

The measurement of the effective Pomeron trajectory at large |t| yields a slope of a' = — 0.0135± 
0. 007 A(stat.) ± 0.0051 (syst.) GeV -2 . This is lower than that observed for elastic J/ip pho- 
toproduction at low \t\ POl and also lower than the slope obtained from hadronic scattering 
{a' = 0.26 ± 0.02 GeV -2 iBTI). The observation of the effective slope being small is compati- 
ble with the predictions of models based on BFKL evolution ifTTl . 

The spin density matrix elements of the J/if; have been extracted in three regions of t. The 
results are found to be consistent with s-channel helicity conservation within the experimen- 
tal uncertainties and, therefore, are compatible with models II11I1IT11 172111311131 in which the 
longitudinal momentum of the photon is shared symmetrically between the quarks of the J/i/j . 
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\t \ range 
(GeV 2 ) 


(\t\) 

\r\/ 

(GeV 2 ) 


da/dt 
(nb/GeV 2 ) 


9 3 
z — o 


9 A'? 


o.iu m u.zy m u.uo 


o — ^ 


^ 4^ 


q ns -i- n 9^ -i- n 


— o 




1 47 + 01^ + 018 


o — u 




u.o i m u.iz m u.ii 


— 7 


(S 47 


n fil n 4- D DQQ 4- D D74 


7-9 


7.92 


0.285 ± 0.046 ± 0.034 


9-12 


10.4 


0.151 ±0.026 ±0.017 


12 - 15 


13.4 


0.093 ±0.020 ±0.010 


15 - 21 


17.7 


0.0236 ±0.0067 ±0.0027 


21-30 


25.0 


0.0045 ± 0.0023 ± 0.0005 



Table 1: The differential photoproduction cross section da/dt in the kinematic range 
50 < W 1P < 150 GeV and z > 0.95. The first uncertainty is statistical and the second is 
systematic. 



W yp range 


(W„) 




(GeV) 


(GeV) 


(nb) 


50-68 


58.4 


7.26 ±0.57 ±0.85 


68-86 


76.5 


8.11 ±0.68 ±0.90 


86 - 104 


94.6 


9.22 ± 0.87 ± 1.06 


104 - 122 


113 


13.5 ± 1.4 ± 1.7 


122 - 140 


131 


13.0 ± 1.8 ± 1.9 


140 - 160 


150 


14.0 ±2.2 ±2.4 



Table 2: The photoproduction cross section as a function of W JP integrated over the kinematic 
range 2 < \t\ < 5 GeV 2 and z > 0.95. The first uncertainty is statistical and the second is 
systematic. 



W JP range 


(W 1P ) 




(GeV) 


(GeV) 


(nb) 


50 - 82.5 


64.4 


1.24 ±0.18 ±0.14 


82.5 - 115 


97.4 


2.75 ±0.35 ±0.31 


115 - 147.5 


130 


3.98 ±0.69 ±0.57 


147.5 - 180 


163 


3.26 ±0.98 ±0.58 



Table 3: The photoproduction cross section as a function of W 1P integrated over the kinematic 
range 5 < \t\ < 10 GeV 2 and z > 0.95. The first uncertainty is statistical and the second is 
systematic. 
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W 1P range 


(W JP ) 




(GeV) 


(GeV) 


(nb) 


50 - 100 


71.0 


0.499 ± 0.093 ± 0.060 


100 - 150 


122 


0.94 ±0.19 ±0.13 


150 - 200 


173 


1.62 ±0.52 ±0.38 



Table 4: The photoproduction cross section as a function of W JP integrated over the kinematic 
range 10 < \t\ < 30 GeV 2 and z > 0.95. The first uncertainty is statistical and the second is 
systematic. 



\t\ range 
(GeV 2 ) 


(1*1) 
( GeV 2 ) 


5 


a 


2-5 


3.06 


0.77 ±0.14 ±0.10 


1.193 ±0.035 ±0.025 


5-10 


6.93 


1.29 ±0.23 ±0.16 


1.322 ±0.057 ±0.040 


10-30 


16.5 


1.28 ±0.39 ±0.36 


1.322 ±0.097 ±0.090 



Table 5: The value of 5 obtained when applying a fit to the data of the form a{W 1P ) oc W 1P for 
each \t\ range, together with the corresponding value of a obtained from a — (S + 4)/4. The 
first uncertainty is statistical and the second is systematic. 



<l*l> 

(GeV 2 ) 


r 04 
' 1-1 


r 04 
'00 


Re{r? 4 } 


3.06 


-0.047 ±0.067 ±0.009 


0.01 ±0.12 ±0.04 


0.022 ±0.069 ±0.035 


6.93 


-0.07 ±0.14 ±0.07 


-0.03 ±0.17 ±0.02 


0.06 ±0.12 ±0.05 


16.5 


-0.19 ±0.22 ±0.12 


0.04 ±0.28 ±0.04 


-0.08 ±0.19 ±0.08 



Table 6: The spin density matrix elements for the kinematic range 50 < W JP < 150 GeV and 
z > 0.95. The first uncertainty is statistical and the second is systematic. The data are quoted 
at the average |t| values over the ranges given in table 
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Figure 1: The invariant mass distribution in the kinematic region 50 < W 1P < 150 GeV, 

z > 0.95 and \t\ > 2 GeV 2 . The histogram shows the sum of the Monte Carlo simulations of 
J / i/j production using HITVM (open histogram), the contribution from lepton pair production as 
simulated by the LPAIR program (dark shaded histogram) and the contribution from diffractive 
ip(2S) events as simulated with the DIFFVM program (light shaded histogram). 
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Figure 2: Kinematic distributions of the dimuon sample in the mass range 2.9 < M M + M - < 
3.3 GeV. a) The polar angle 9^ and b) the transverse momentum p% of the muon tracks, c) 
The elasticity z and d) the photon-proton centre-of-mass energy W 1V . e) The distribution of the 
cosine of the polar angle and f) the azimuthal distribution of the positively charged decay muon 
in the helicity frame, g) The distribution of the squared dimuon transverse momentum. 
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Figure 3: The photon-proton differential cross section dcr/dt for J/tp production in the kine- 
matic range 50 < W w < 150 GeV, \t\ > 2.0 GeV 2 and z > 0.95. The inner error bars 
correspond to the statistical error and the outer error bars are the statistical and systematic er- 
rors added in quadrature. The solid line shows the prediction from the BFKL calculation in the 
leading logarithmic approximation for fixed a s lfT51 . The dashed (dotted) curve corresponds to 
the BFKL calculation including non-leading corrections and using a fixed (running) a s lfT5ll . 
The dashed-dotted curve, shown in the range \t\ < M%^, shows a calculation based on the 
DGLAP equation in the leading logarithmic approximation lfT3l . 
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Figure 4: The photon-proton cross section as a function of W JP in three bins of \t\. The in- 
ner error bars correspond to the statistical error and the outer error bars are the statistical and 
systematic errors added in quadrature. The solid lines show the predictions from the BFKL 
calculation in the leading logarithmic approximation and the dashed lines correspond to the 
BFKL calculation including non-leading corrections using a fixed a s ITT31 . The dashed-dotted 
curve is the result of a calculation based on the DGLAP equation in the leading logarithmic 
approximation lfl3l . 
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Figure 5: Normalised decay angular distributions for J ftp meson production in three bins of 
|f|: a,b) 2 < \t\ < 5 GeV 2 ; c,d) 5 < \t\ < 10 GeV 2 and e,f) 10 < \t\ < 30 GeV 2 . The left 
column (a,c,e) shows the azimuthal distributions of the positively charged decay muon in the 
helicity frame and the right column (b,d,f) shows the polar angle distributions. The inner error 
bars show the statistical error and the outer error bars show the statistical and systematic errors 
added in quadrature. The solid lines show the results of a two-dimensional fit to the data (see 
text). The shaded band represents the statistical uncertainty for the fit. The dashed line shows 
the expectation from s-channel helicity conservation. 
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Figure 6: The three spin density matrix elements a) r^_ x , b) Tqq and c) Re{r°o} for the J/ip as a 
function of The inner error bars represent the statistical uncertainty and the outer error bars 
the combined statistical and systematic uncertainties. The dashed line shows the expectation 
from SCHC. The results from the ZEUS collaboration for the photoproduction of J/tp and p° 
mesons ifHI are also shown. 
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